Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 5 February 2008 (MN WF$t style file v2.2) 



Intermediate-mass black holes in dwarf galaxies: the case 
of Holmberg II 
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ABSTRACT 

In order to constrain the density of intermediate-mass black holes (IMBHs) in galaxies, 
we run smoothed particle hydrodynamics (SPH) simulations of a gas-rich disc dwarf 
galaxy, where different halo and disc populations of IMBHs are embedded. IMBHs, 
when passing through dense gas regions, can accrete gas and switch on as X-ray 
sources. We derive the luminosity distribution of simulated IMBHs, by assuming that 
they accrete at the Bondi-Hoyle rate. The X-ray distribution of simulated IMBHs 
has been compared with that of observed sources in the dwarf galaxy Holmberg II, 
chosen for its richness in gas, its small mass (compared to spiral galaxies), and the 
accuracy of the available X-ray measurements. Holmberg II also hosts one of the 
strongest IMBH candidates. From this comparison, we find that the density parameter 
of disc (halo) IMBHs must be ft, £ 10~ 5 Oj, (fi. ^ 10~ 2 £lfc, where il b is the density 
parameter of baryons), for a radiative efficiency 10 -3 and an IMBH mass of 10 4 M Q . 
These constraints imply that a dwarf galaxy like Holmberg II cannot host more than 
1 (1000) disc (halo) 10 4 M Q IMBHs. 

Key words: black hole physics - methods: A^-body simulations - galaxies: individual: 
Holmberg II - X-rays: general 



1 INTRODUCTION 

The existence of intermediate-mass black holes (IMBHs), i.e. 
black holes (BHs) with mass from ~ 20 to 1O 5 M , is still 
controversial (see van der Marel 2004 for a review). Various 
IMBH formation processes have been proposed: (i) collapse 
of massive ( £ 260 M®) metal free stars (Heger & Woosley 
2002); (ii) runaway collision of stars in young clusters (Porte- 
gies Zwart & McMillan 2002); (iii) repeated mergers of stel- 
lar mass BHs in globular clusters (Miller & Hamilton 2002) , 
and (iv) direct collapse of dense, low angular momentum gas 
(Haehnelt & Rees 1993; Umemura, Loeb & Turner 1993; 
Loeb & Rasio 1994; Eisenstein & Loeb 1995; Bromm & 
Loeb 2003), driven by turbulence (Eisenstein & Loeb 1995) 
or gravitational instabilities (Koushiappas, Bullock & Dekel 
2004; Begelman, Volonteri & Rees 2006; Lodato & Natara- 
jan 2006). 

Recent spectroscopic and photometric measurements 
suggest the presence of IMBHs in globular clusters (Geb- 
hardt, Rich & Ho 2002, 2005; Gerssen et al. 2002; van den 
Bosch et al. 2006). Furthermore, IMBHs accreting dense 
gas (Krolik 2004; Mii & Totani 2005) or matter from com- 
panion stars (Hopman, Portegies Zwart & Alexander 2004; 
Kalogera et al. 2004; Portegies Zwart, Dewi & Maccarone 



2004; Hopman & Portegies Zwart 2005; Patruno et al. 2005; 
Baumgardt et al. 2006) have been invoked to explain the 
presence of ultra- luminous X-ray sources (ULXs), i.e. X-ray 
sources with luminosity higher than 10 39 erg s _1 (see Col- 
bert & Miller 2005 for a review). 

In general, IMBHs accreting gas in molecular clouds 
are expected to switch on as X-ray sources (either ultra- 
luminous or not). Mapelli, Ferrara & Rea (2006; hereafter 
MFR) calculated, via dedicated N-body simulations, the 
number of X-ray sources which are expected to originate 
from IMBHs accreting gas in the Milky Way. They com- 
pared the number of simulated sources with the number of 
unidentified X-ray sources observed in our Galaxy, in order 
to find an upper limit on the number of IMBHs. The results 
of MFR depend on the adopted accretion model and on the 
spatial IMBH distribution. In a conservative case, i.e. assum- 
ing that the accretion mechanism on IMBHs is an advection 
dominated accretion flow (ADAF) with radiative efficiency 
r] = 10~ 3 , the density of 10 4 M© IMBHs in our Galaxy must 
be £ 10 -1 fif, (where = 0.042 is the density of baryons in 
terms of the critical density of the Universe, Spergel et al. 
2006), if the IMBHs are distributed according to a Navarro 
Frenk & White (1996; NFW) profile. Diemand, Madau & 
Moore (2005; DMM) recently found that objects formed in 
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high u-fluctuations follow a more concentrated distribution. 
Adopting the DMM profile, the upper limit on the IMBH 
density becomes ^ 10~ 2 f2f,. 

In this paper we adopt the same technique developed 
in MFR (with some substantial improvements, such as the 
treatment of the gas and of the dark matter halo), in order 
to put constraints on the IMBH population of Holmberg 
II (=UGC 4305=DDO 50; hereafter HoII), a dwarf irregu- 
lar galaxy belonging to the M81 group (Karachentsev et al. 
2002). We chose a dwarf galaxy, and in particular HoII, for 
various reasons. First of all, there are many papers aimed to 
constrain the number of IMBHs in the Milky Way, from a dy- 
namical point of view (Carr & Sakellariadou 1999; see also 
Lacey & Ostriker 1985; Wasserman & Salpeter 1994; Mu- 
rali, Arras & Wasserman 2000 and reference here), from the 
comparison with X-ray sources (Mii & Totani 2005; MFR) 
and from the study of globular cluster formation and disrup- 
tion (Ostriker, Binney & Saha 1989; Moore 1993; Klessen & 
Burkert 1996; Arras & Wasserman 1999; Murali et al. 2000); 
whereas there are no papers dedicated to constrain IMBHs 
in dwarf galaxies. Is the density of IMBHs comparable in 
dwarf and normal galaxies? A dependency of the density of 
IMBHs on the mass of the host galaxy could have important 
implications on the theoretical models of IMBH formation. 

Furthermore, dwarf galaxies are much less massive than 
the Milky Way, allowing us to reach higher mass and spatial 
resolution in our simulations (~ 10 2 Mq and ,$ 1 pc, re- 
spectively, in the case of HoII, versus 5 x 10 4 Mq and ~ 100 
pc for the Milky Way; see MFR). 

Among the possible candidates, HoII is particularly 
suitable for our analysis, because it is a very gas-rich dwarf 
galaxy, about one tenth of its kinematic mass being repre- 
sented by HI (Puche et al. 1992; Stewart et al. 2000; Bu- 
reau & Carignan 2002). Therefore, the probability that an 
IMBH passes through a dense gas region is much higher in 
HoII than in the Milky Way. In addition, the rotation curve 
(and, consequently, the dynamical mass) and the other dy- 
namical properties of HoII have been accurately measured 
in the last years (Bureau & Carignan 2002 and references 
there). Finally, the X-ray sources detected in HoII have been 
reasonably well studied (Kerp, Walter & Brinks 2002), espe- 
cially X-l, which is one of the most powerful observed ULXs 
(Dewangan et al. 2004) and one of the strongest IMBH can- 
didates. By the way, the presence of this ULX ensures that, 
if IMBHs exist and are connected with ULXs, HoII is one of 
the best places where we can find them. 

More details about the observed properties of HoII are 
given in Section 2. In Section 3 we describe our simulations. 
The results have been presented and compared with observa- 
tions in Section 4 (for IMBHs accreting molecular hydrogen) 
and Section 5 (atomic hydrogen). Section 6 is a summary of 
the main findings. 



2 OBSERVATIONAL PROPERTIES OF HOII 

HoII is a very gas-rich (Puche et al. 1992; Stewart et al. 2000; 
Bureau & Carignan 2002) dwarf irregular, belonging to the 
M81 group (Karachentsev et al. 2002). Its kinematic total 
mass, derived from HI measurements (Bureau & Carignan 
2002), is 6.3 x 1O 9 M , and it is mainly dark (~ 80 per cent). 
The mass in HI is 6.44 x 10 8 Mq, indicating that there is more 



luminous mass in gas than in stars. The circular velocity, 
~ 36 km s _1 (see Bureau & Carignan 2002 and Fig. 1), 
implies that this galaxy rotates quite slowly. There are no 
signs of interaction between HoII and close-by members of 
the M81 group (Kar 52 and UGC 4483), and HoII is only 
slightly affected by ram pressure (in form of a faint HI tail; 
Bureau & Carignan 2002). So it makes sense to simulate 
HoII as an isolated galaxy. 

Kerp et al. (2002), by using ROSAT data, detect 31 
X-ray sources (with luminosity Lx Si 10 37 erg s _1 ) located 
in HoII. To avoid contamination from background sources, 
they consider only those sources which have been identified 
in at least one additional frequency range. By adopting this 
selection criterion, Kerp et al. (2002) are left with 13 sources, 
which should belong to HoII. Some of them can be associ- 
ated with either X-ray binaries or supernova remnants (Kerp 
et al. 2002), but this identification is not certain. Then, in 
this paper, we will make the very conservative assumption 
that all these 13 sources are unidentified and could be, in 
principle, IMBHs. 

The brightest among these X-ray sources, HoII X-l, is 
one of the most powerful ULXs observed up to now, having 
a X-ray luminosity Lx ~ 2 x 10 40 erg s _1 in the brightest 
state (Dewangan et al. 2004). An ionized nebula has been 
detected around HoII X-l (Pakull & Mirioni 2003; Kaaret, 
Ward & Zezas 2004; Lehmann et al. 2005), and important ra- 
dio emission has been shown to be connected with it (Miller, 
Mushotzky & Neff 2005). The spectrum of HoII X-l can be 
fitted combining an absorbed power law and a thermal com- 
ponent (thermal plasma or multi-colour blackbody), which 
indicates relatively low temperatures (~ 0.14 — 0.22 keV; 
Lehmann et al. 2005). For all these reasons, the hypothesis 
that HoII X-l is a beamed X-ray binary is disfavoured, and 
this source is considered one of the strongest IMBH candi- 
dates (Miller, Fabian & Miller 2004). 



3 NUMERICAL SIMULATIONS 

As in MFR, the simulations have been carried out using the 
parallel N-body code GADGET-2 (Springel 2005). The sim- 
ulations were performed using 8 nodes of the 128-processor 
cluster Avogadro at the Cilea (http:/ /www. cilea.it). We sim- 
ulated a N-body model of HoII, in which we embed a halo 
population of IMBHs. 



3.1 HoII model 

To simulate a model of HoII, we combined three different 
components: a NFW halo, a stellar exponential disc and a 
uniform molecular gas disc (or an exponential atomic gas 
disc). 

Halo, disc and gas particle velocities are generated using 
the Gaussian approximation (Hernquist 1993), as described 
in MFR. In the following, we describe the details of our 
model for halo, stellar disc and gas. All the main adopted 
parameters of this model are listed in Table 1. 
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corresponding mass of dark matter within -R200 is Mdm 




Figure 1. Rotation curve of HoII. The solid (dashed) line refers 
to the average circular velocity of disc star particles assuming a 
halo concentration c = 3 and a uniform disc of molecular gas (ex- 
ponential disc of atomic hydrogen). Data points are from Bureau 
& Carignan (2002), and refer to HI measurements. 



3.1.1 Halo 

The parameters of the NFW halo have been derived as de- 
scribed in MFR. In particular, the density distribution of 
the halo is (NFW; Moore et al. 1999): 

ps 



Ph( - r) {r/r s )~t\l + (r/r a )°']V-1)/<* , 

where we choose (a,/3, 7) = (1,3, 1), and p s — p C ritS c , p, 
being the critical density of the Universe and 

200 



(1) 



(2) 



3 lrr(l+c)-[c/(l+c)]' 

where c is the concentration parameter and r s is the halo 
scale radius, defined by r s = R200/C R200 is the radius en- 
compassing a mean overdensity of 200 with respect to the 
background density of the Universe. -R200 can be calculated 
as -R200 = V2oo/[10 H(zvir)], where V200 is the circular ve- 
locity at the virial radius and H{z vir ) is the Hubble param- 
eter at the virialization redshift z V i r . For HoII, we adopt 
V200 = 36 km s _1 (Bureau & Carignan 2002) and z V i r = 1.8, 
obtaining -R200 — 20 kpc, in agreement with HI observations 
(Bureau & Carignan 2002). In order to reproduce the ob- 
served rotation curve (Fig. 1), we have to assume c < 3, 
which is a very low concentration. Cosmological simulations 
(Bullock et al. 2001) show that dwarf galaxies should have 



significantly higher concentration 



10). However such 



high concentrations cannot account for the observed rota- 
tion curve, indicating that the NFW profile fails to predict 
the properties of dwarf irregular galaxies. For the problem of 
fitting dwarf galaxy rotation curves by using a NFW profile 
see Navarro (1998) and Mo & Mao (2000). 

Adopting the parameters listed above, the total mass 
within -R200 is M200 = 5.9 x 1O 9 M0, in nice agreement with 
observations (6.3 x 10 9 M Q , Bureau & Carignan 2002). The 



M 2 , 



M d - M 



where Md is the stellar disc mass, and 



M g as is the gas mass (M gaa is equal to the mass in molecular 
gas, M mo i, or to the mass in atomic gas, Mm, depending on 
whether we simulate molecular or atomic gas, see Section 
3.1.3). 



3.1.2 Disc 

The stellar disc profile is (Hernquist 1993): 

M d 



pd(R,z) = 



-R/R-d sech 2 (z/zo), 



(3) 



where Md — 4 x 10 s Mq is the disc mass, R d = 2 kpc the 
disc scale length and zq = 0.1 Rd the disc scale height. These 
adopted parameters are in good agreement with observa- 
tions (Bureau & Carignan 2002) and correspond to a spin 
parameter A ~ 0.07, indicating a dynamically stable disc 
(Mo, Mao & White 1998). 

3.1.3 Gas component 

The gas in HoII is expected to be partially molecular and 
partially atomic. We are particularly interested in the molec- 
ular gas, because it is expected to be more efficiently ac- 
creted by IMBHs (being much denser than the atomic hy- 
drogen). However, we want to check the possibility that also 
IMBHs accreting atomic hydrogen are observable as X-ray 
sources (see MFR for the Milky Way). To simulate these 
two gas components simultaneously is beyond the purpose 
of this paper. Thus, we made two different sets of runs (see 
Section 3.3), part of them including only molecular gas and 
the remaining only atomic gas. In runs with molecular gas 
the mass of gas which is expected to be in atomic form is 
attributed to dark matter particles. 

We chose the profile of the gas component according 
to whether the gas is assumed to be molecular or atomic. 
The molecular gas has been represented as a uniform disc 
with scale length R mo i = 500 pc, scale height z mo i = 60 pc 
(Taylor et al. 1999) and mass M moi = 2 x 10 s M Q . In fact, 
HoII has a total HI mass 6.44 x 1O 8 M (Bureau & Cari gnan 
2002) and we know that the average fraction of H2 versus 
HI in dwarf irregular galaxies is ~ 0.3 (Leroy et al. 2005). 

The adopted equation of state is polytropic: 



Rgas — K Pgas 



(4) 



where P gas and p ga s are the gas pressure and density, V is 
equal to 5/3 (which avoids excessive fragmentation, ensur- 
ing the stability of the disc, Escala et al. 2005) and k is 
a parameter corresponding to the entropy of the gas. k is 
related to the gas temperature by the following relation. 



■ 60K 



9 x lO^g-^" 1 ) cmV 



( 



r ( 



V 100cm" 3 



where n gas is the gas density and p the mean molecular 
weight (p — 2 for molecular gas). We consider values of k 
ranging from 9 x 10 23 g -2 ^ 3 cm 4 s -2 (corresponding to a gas 
temperature T gas — 60 K, suitable for molecular gas; see 
Spitzer 1998) to 1.5 x 10 25 g" 2/3 cm 4 s" 2 (corresponding to a 
gas temperature T ga3 — 1000 K, as an upper limit). Discs 
with a low k tend to become clumpy before hotter discs 
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Table 1. Initial parameters for the HoII model. 
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v 2 


36 km s~ 1 


R200 


19.6 kpc 


M200 


5.9xlO 9 M 


A 


0.07 


M d 


4xl0 8 M o 


Rd 


2 kpc 


20 


0.1 R d 


M mol 


2xlO 8 M 




500 pc 


Z"tnol 


60 pc 


M m 


6.44xl0 8 M o 


Rm 


5 kpc 


ZHI 


600 pc 



(Escala et al. 2005) , and have a smaller Jeans mass. On the 
other hand, our choice of F = 5/3 is very conservative, and 
balances the effect of k, avoiding extreme fragmentation of 
the disc. 

Instead, when we simulate atomic gas, we assume an 
exponential distribution (the same as equation (3)), accord- 
ing to the observed neutral hydrogen exponential profile in 
disc galaxies (Lockman 2002) . The total mass of atomic hy- 
drogen is assumed to be Mm = 6.44 x 10 8 Mq (Bureau & 
Carignan 2002), the scale length of the disc is Rm = 5 kpc, 
and the scale height zm = 600 pc (Table 1). The equation 
of state assumed for the neutral atomic gas is the same as 
that for the molecular gas (with /1 = 1.3, T gas = 60 K and 
rigas = 1 cm -3 ). 

In the simulations with atomic exponential gas, we 
adopt the same values of c, -R200, M200, Md, Rd and zo (Ta- 
ble 1) that we used for the runs with molecular gas. In fact, 
the rotation curve is nearly the same for runs with molec- 
ular and with atomic hydrogen (Fig. 1), because the gas 
component (about one tenth of the dark matter mass) has 
a minor impact on the dynamics. Thus, the only differences 
between runs with molecular and with atomic hydrogen are 
the treatment of the gas component and the total mass in 
dark matter 1 . This choice makes easier the comparison be- 
tween the two cases, since the main characteristics of the 
galaxy model are unchanged. 

3.2 Intermediate-mass black holes 

The number of IMBHs hosted in galaxies is completely un- 
known. If we assume that they formed at high redshift, then 
an indicative estimate of the IMBH number has been de- 
rived in Volonteri, Haardt & Madau (2003). As in MFR, 
we use the formalism of Volonteri et al. (2003) to generate 
the initial conditions of our simulations. In particular, we 
estimate the number of IMBHs residing in a galaxy (N.) as 

< 6 > 

1 When we run simulations with molecular gas but without 
atomic hydrogen, we have to redistribute among dark matter par- 
ticles the mass which is expected to be in form of atomic hydrogen 
(see Section 3.3). 



where fi. and Qb (=0.042; Spergel et al. 2006) are the current 
density of IMBHs and baryons, respectively, in terms of the 
critical density of the Universe; M b is the mass in baryons 
of the host galaxy, and m. is the current average IMBH 
mass. For HoII Mb ~ 10 9 Mq, so that equation (6) becomes 
N. ~ 100 (n./10" 3 n 6 ) (10 4 M e /m.). According to MFR, 
the case with Q,. — 10~ 3 f2(, corresponds to assume that 
IMBHs formed in 3 a fluctuations at redshift z ~ 20 — 25. 

This estimate is affected by a number of uncertainties, 
and we adopt it only because it is a straightforward formal- 
ism to express the number of IMBHs in terms of the critical 
density of the Universe. However, the results of our simu- 
lations have more general validity, and can be applied also 
to IMBHs formed via different mechanisms (e.g. runaway 
collapse) . 

The spatial and velocity distribution of IMBHs are an- 
other severe problem. If IMBHs formed in high redshift mini- 
halos, as remnants of metal free stars (Volonteri et al. 2003), 
they probably constitute a halo population in current galax- 
ies. For this reason, Mii & Totani (2005) chose a NFW pro- 
file to represent the distribution of IMBHs in their model. 
However, White & Springel (2000) suggested that remnants 
of population III stars should be more concentrated inside 
present-day halos than younger objects. DMM confirmed 
this idea and showed that the current spatial distribution of 
objects formed in high-cr fluctuations depends only on the 
rarity of the peak in which they are born. DMM rewrote the 
NFW profile, accounting for this correction, as 



(r/r„)T [1 + (r/r„)°]tf"--r)A« ' V ; 

where a and 7 are the same as defined in the previous sec- 
tion; r„ = r 3 /f v is the scale radius for objects formed in a 
v a fluctuation (with f v = exp (u/2)), and /3„ = 3+0.26 v x % . 
Hereafter, we will refer to equation (7) as the DMM profile. 
This profile is likely to track the distribution of IMBHs, if 
they are a halo population (MFR). 

On the other hand, IMBHs could also be a disc pop- 
ulation, either if they form by runaway collapse in young 
clusters, or if they are remnants of population III stars born 
in the disc. The latter scenario is not so unrealistic, if a part 
of the metal free gas (i.e. Z ^ 10~ 4 , Bromm et al. 2001; 
Schneider et al. 2002) remains unpolluted even at low red- 
shift, forming population III stars in the disc of galaxies 
(Krolik 2004). For example, Jimenez & Haiman (2006) ar- 
gued that a percentage of metal free stars as high as 10-50 
per cent of the total stellar population is required at redshift 
z ~ 3 — 4 to explain some puzzling observations (e.g. the sig- 
nificant ultraviolet emission from Lyman Break Galaxies at 
wavelengths shorter than 912 A). If IMBHs are a disc pop- 
ulation, they are so massive that dynamical friction exerted 
by gas is expected to pull them toward the galaxy center in 
a short time-scale (~3x 10 s yr, for 10 4 M© IMBHs, if the 
average density of gas in the stellar disc is ~ 1 cm -3 ; Binney 
& Tremaine 1987). Then, we expect that their distribution 
is more concentrated than that of other disc populations. 

We will examine both these cases, i.e. IMBHs dis- 
tributed in the halo and in the disc. For the halo IMBHs, we 
adopt a DMM profile, as described in MFR. In MFR we also 
discussed a NFW distribution; but this distribution is quite 
unlikely for halo IMBHs (White & Springel 2000; DMM). 
Then, we decided not to consider this case in the present 



IMBHs in dwarf galaxies 5 



paper. Furthermore, MFR showed that the constraints for 
NFW IMBHs have the same trend of those derived for DMM 
IMBHs, but are a factor of 10 weaker. For disc IMBHs, we 
assume an exponential distribution, similar to that described 
for stars, but with a cut-off radius R c = 1 kpc, a scale ra- 
dius Rd, bh = 500 pc and a scale height Rd, bh = 10 pc, 
accounting for dynamical friction. 

3.3 Description of runs 

The simulated dwarf galaxy is composed by: 

- 10 6 halo particles inside -R200, each one of 5.3 x 10 3 M Q 
or 4.9 x 10 3 M (in the case we simulate a uniform disc of 
H2 or an exponential disc of atomic gas, respectively); 

- 10 disc particles, each one of 4 x 10 2 M Q ; 

- 10 6 smoothed particle hydrodynamics (SPH) gas parti- 
cles, each one of 2 x 10 2 M or 6.4 x 10 2 M Q (in the case 
we simulate a uniform disc of H2 or an exponential disc of 
atomic gas, respectively). 

The gravitational softening lengths are l so ft = 6, 0.2 and 
0.1 pc for halo, disc and gas particles, respectively. The av- 
erage initial SPH length 2 is ?sph ~ 2 pc. 

It is worth noting that (i) we simulated a 'living' halo, 
rather than the rigid halo adopted by MFR; (ii) we make an 
SPH treatment of gas particles, whereas MFR used a Monte 
Carlo technique to account for molecular and atomic gas in 
their merely N-body simulations. 

In this galaxy model we followed the dynamics of vari- 
ous populations of IMBHs. The mass of each IMBH is 10 4 or 
10 5 Mq depending on the run. We also made some checks 
with higher mass BHs (10 6 ~ 7 M Q ). We did not consider 
lower mass BHs, because of our mass resolution. The num- 
ber of IMBHs per simulation is 10—1000 (corresponding to 
n. = 10~ 4 — 10~ 2 fi 6 ). IMBHs were distributed according 
to a DMM profile or an exponential disc. The softening of 
IMBH particles is as small as possible (0.002 pc). 

We made more than 30 different runs; but, for sim- 
plicity, in Table 2 we list only the most significant among 
them. In eight of the runs reported in Table 2 (disc4_10, 
disc4_100, disc5_10, K_disc4_100, DMM4.100, DMM4.1000, 
DMM5.100 and DMM5_1000) we simulated a uniform disc 
of molecular gas; while in the remaining three (HI_disc4_10, 
HI_DMM4_100 and HI_DMM5_100, i.e. the runs labelled as 
'HI') we simulated an exponential disc of atomic hydrogen. 

IMBHs are assumed to be a disc population in five 
of these runs (disc4_10, disc4_100, disc5_10, K_disc4_100 
and HI_disc4_10, i.e. the runs labelled as 'disc') and a 
halo population, following a DMM profile, in the other six 
runs (DMM4.100, DMM4.1000, DMM5.100, DMM5.1000, 
HI_DMM4_100 and HI_DMM5_100, i.e. the runs labelled as 
'DMM'). Finally, those runs whose label ends in 4_n (5_n) 

2 Bate & Burkert (1997) showed that, if Isph ~ lj, where lj is 
the Jeans length, and if igpn S> I soft or ^SPH l so jt, some un- 
physical effects can be induced. In our simulations we always have 
lj ~ 10 — 30 pc > Zsph- Thus, even if the balance between pres- 
sure and gravity at distances smaller than ~ 2 pc can be affected 
by our choice of isPH and l so ft, the Jeans mass is resolved by 
the simulations. Furthermore, we compared our simulations with 
some test runs with l ao ft = 6, 2 and 2 pc (=£sph) for halo, disc 
and gas particles, respectively, and we found negligible differences. 



host a number n (with n=10, 100 or 1000) of 10 4 M (10 5 
M ) IMBHs. 

The gas has n = 9 x 10 23 g~ 2/3 cm 4 s -2 (see equa- 
tion 4) in all the cases reported in Table 2, except for the 
run K_disc4_100, where k = 1.5 x 10 25 g~ 2/3 cm 4 s~ 2 . 

Each run was stopped after 6 Myr. A longer time would 
require an excessive amount of computational time (i.e. more 
than 2 weeks on 16 3-GHz CPUs per run). However, 6 Myr 
is a sufficiently long time, if compared to the expected duty 
cycle of such sources. In fact, IMBHs are expected to rapidly 
heat the surrounding gas and enhance the pressure, dramat- 
ically quenching the accretion rate on a time scale shorter 
than 6 x 10 5 yr (Krolik 2004). 



4 IMBHS ACCRETING MOLECULAR GAS 

The aim of this paper is to constrain the number of IMBHs 
in HoII following the same technique adopted for the Milky 
Way by MFR. The basic idea is that, if a certain number 
of IMBHs are hosted in HoII and pass through dense gas 
regions, they should accrete this gas and possibly become 
observable as X-ray sources. 

Unfortunately, we do not know the details of the accre- 
tion process for IMBHs. As the simplest approximation (Mii 
& Totani 2005; MFR), we can assume that these IMBHs ac- 
crete at the Bondi-Hoyle luminosity: 

L x {p g , v) = Anrjc 2 G 2 m, 2 p g v~ 3 , (8) 

where r\ is the radiative efficiency, c the speed of light, G the 
gravitational constant, m. the IMBH mass, p g the density 
of the gas surrounding the IMBH. v = (v 2 + o 2 MC + c 2 ) 1//2 , 
where v is the relative velocity between the IMBH and the 
gas particles, ctmc and c s are the molecular cloud turbu- 
lent velocity and gas sound speed, respectively. We adopt 
&mc = 3.7 km s _1 (MFR and references therein), and we 
extract v and c s directly from our simulations (see next sec- 
tion). Lx in the above equation indicates the total X-ray lu- 
minosity. Properly speaking, the Bondi-Hoyle formula refers 
to the bolometric luminosity. However, in the case of ULXs 
the X-ray luminosity is much higher than the optical lumi- 
nosity (Winter, Mushotzky & Reynolds 2006), justifying our 
approximation. 

The most slippery quantity in equation (8) is the effi- 
ciency 77, which strongly depends on the details of the ac- 
cretion process. Agol & Kamionkowski (2002) show that 
density gradients in the molecular cloud are likely to im- 
print a small amount of angular momentum on the accret- 
ing gas. This allows the formation of an accretion disc for 
most of IMBHs. However, this disc is not necessarily thin. 
It is more likely that IMBHs are surrounded by a thick disc, 
e.g. an ADAF disc (Narayan, Mahadevan & Quataert 1998; 
Quataert & Narayan 1999). The radiative efficiency of thick 
discs is quite lower than for thin discs. 

In particular, the ADAF luminosity is oc M 2 (where M 
is the accretion rate) instead of oc M as for the Bondi-Hoyle 
luminosity. However, for the systems we are considering 




where M Edd is the Eddington accretion rate and n gaa the 
gas density. For the range M/Msdd ~ 10 _1 — 10~ 2 it is 
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Table 2. Initial parameters for IMBHs and gas distribution. 



Run IMBH mass (M ) Number of IMBHs H./Qf, IMBH profile Gas profile 



disc4_10 


10 4 


10 


io- 


4 


exponential disc 


uniform disc 


disc4_100 


10 4 


100 


io- 


3 


exponential disc 


uniform disc 


disc5.10 


10 5 


10 


io- 


3 


exponential disc 


uniform disc 


K_disc4_100 


10 4 


100 


io- 


3 


exponential disc 


uniform disc 


DMM4.100 


10 4 


100 


10" 


3 


DMM 


uniform disc 


DMM4.1000 


10 4 


1000 


10" 


2 


DMM 


uniform disc 


DMM5.100 


10 5 


100 


io- 


2 


DMM 


uniform disc 


DMM5.1000 


10 5 


1000 


io- 


1 


DMM 


uniform disc 


HI_disc4_10 


10 4 


10 


io- 


4 


exponential disc 


exponential disc 


HI_DMM4_100 


10 4 


100 


io- 


3 


DMM 


exponential disc 


HI.DMM5.100 


10 5 


100 


10- 


2 


DMM 


exponential disc 



Figure 2. Face-on projection of the central 200 pc of the simulated gas disc after 3 Myr. The color coding shows the z-averaged gas 
density. The scale, which is logarithmic, approximately goes from 0.1 (black; blue in the online version) to 10 4 (white; bright yellow in 
the online version) cm -3 . The bright yellow spots correspond to high-density regions hosting IMBHs. Left panel: run disc4_100. Right 
panel: DMM4.100. The number of accreting IMBHs is clearly higher for the case disc4_100 than for DMM4.100. 



reasonable to assume, as a lower limit, r\ — 10 3 (Narayan 
et al. 1998). 

Then, in this paper we will assume tj — 10 -3 as a conser- 
vative reference value. According to equation (8), our results 
can be easily rescaled for different values of rj. 



4.1 Gas density and velocity 

The three terms in equation (8) that we can directly derive 
from our simulations are p g , v and c s . 

The high resolution of our simulations allows us to com- 
pute the gas density p g close to the IMBH accretion radius 



Gm, 

r a cc = — 5— ~ 0.03 pc 



10 4 M F 



( v - 

V40 km s" 1 



(10) 



In particular, in order to derive the density p g required 
in equation (8) , we count the gas particles which fall within 
a certain neighborhood of each IMBH. The radius of the 
neighborhood has been chosen to be r g = 5 pc for molecular 



hydrogen and r g = 100 pc for atomic hydrogen. In fact, 
the adopted softening allows us to assume radii as small 
as 0.1 — 1 pc; but we require each neighborhood to host 
a number of gas particles sufficient ( ^ 10) to smear out 
numerical fluctuations. Then, p g has been derived as 



Pa = 



No 



(11) 



where N 9 is the number of gas particles in each neighborhood 
of radius r g and m g is the mass of a gas particle (m g — 
2 x 10 2 Mq for molecular gas; see Section 3.3). 

Because of the adopted polytropic equation (4), the gas 
is allowed to become clumpy (Escala et al. 2005). There- 
fore, IMBHs particles attract nearby gas particles, rapidly 
increasing their surrounding density. This is particularly ev- 
ident for disc BHs and for a uniform disc of gas. In Fig. 2, the 
two runs disc4_100 (left panel) and DMM4_100 (right panel) 
are compared. The two panels of this figure show the face-on 
projected density of the gas in the central region of the sim- 
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Figure 3. Simulation disc5_10. Number of IMBHs (filled his- 
tograms) and gas particles (open histograms) which are sur- 
rounded by a number of gas particles N 3 (within a sphere of 
radius r g = 5 pc). The gas particles represented here are a sam- 
ple of 100 particles randomly selected from the simulation. Left 
panel: initial conditions. Right panel: after 6 Myr. After 6 Myr the 
number of gas particles N 9 surrounding the IMBHs has increased 
by a factor of ~ 100. 



ulated galaxy after 3 Myr. The bright spots represent the 



highest density clumps 



1(T 



in the simulation. 



Each of these spots hosts a central IMBH. The density of gas 
clumps surrounding the IMBHs in the left panel (where the 
IMBHs are simulated as a disc population) is significantly 
higher than in the right panel (where the IMBHs follow a 
DMM profile). 

Because IMBHs attract nearby gas particles, the uni- 
form gas disc becomes more and more clumpy and the den- 
sity p g in the neighborhoods of the IMBHs increases (see 
Fig. 3), depending on the IMBH mass, position and veloc- 
ity. This has a significant effect on the estimated luminosity 
(see next section for details). 

Since we sampled p g in neighborhoods of radius r g = 5 
pc (r g = 100 pc for atomic gas) surrounding the IMBHs, 
a natural choice for deriving the relative velocity v is to 
calculate the average relative velocity between the IMBH 
and the gas particles hosted in the same neighborhood. 

Similarly, the sound speed of the gas, c s , around each 
IMBH is calculated as the average sound speed of gas par- 
ticles within the same neighborhood of radius r g , by using 
the relation <? a — 2u g (where u g is the average internal en- 
ergy per unit mass of the gas particles within the consid- 
ered neighborhood). The values of c s calculated with this 
method strongly depend on the density within the neigh- 
borhood, thus reflecting the feedback of the IMBH on the 
gas particles. In fact, the IMBH increases the gas density in 
its neighborhood and consequently heats the gas (because 
of the polytropic equation of state), quenching the accretion 
rate. For example, for disc IMBHs (the more efficient in at- 
tracting gas) we find c B ~ 1 — 10 km s _1 (depending on p g ), 



quite higher than the value c s = 0.3 km s _1 estimated for 
unperturbed Galactic molecular clouds (see MFR and ref- 
erences therein). Of course, this is not a realistic treatment 
of the accretion process; but it allows us to approximately 
account for the feedback from the IMBH. 



4.2 Comparison with observed X-ray sources 

We can now estimate the X-ray luminosity of the simulated 
IMBHs directly from these values of p g , v and c s . 

In Fig. 4 we show the cumulative distribution of X-ray 
sources as a function of Lx for the most significant cases 
and assuming r\ = 10~ 3 . The luminosity distributions can 
be easily rescaled for different values of r\ according to equa- 
tion (8). The left panels show the initial conditions, while 
the right panels represent the evolution after 6 Myr. 

The IMBHs have a significantly different behaviour de- 
pending on their distribution (exponential or DMM). Disc 
IMBHs tend to attract gas particles and to increase their 
luminosity by a factor of ~ 2-10 in 6 Myr, depending on 
their mass. 

On the contrary, most of the halo IMBHs do not attract 
efficiently gas particles. Halo IMBHs do not increase the 
density p g and maintain nearly the same luminosity during 
all the simulation. 

This is probably due to the different velocity distribu- 
tions. In fact, disc IMBHs are substantially corotating with 
the gas particles; whereas halo IMBHs have (on average) 
higher relative velocities and periodically leave the gaseous 
disc. 

This result strongly affects any attempt to constrain 
the number of IMBHs. In fact, if IMBHs are a disc popu- 
lation (as predicted, for example, by the runaway collapse 
scenario; Portegies Zwart & McMillan 2002), the constraints 
on their number are very strong (see also Krolik 2004). 
Even 10 IMBHs of 10 4 M Q , corresponding to Q. = 10~ 4 n b , 
would overproduce the number of X-ray sources in the 
range from 10 38 to 10 40 erg s _1 , unless their radiative effi- 
ciency is 10~ 4 (i.e. a factor of 10 lower than the assumed 
r/^lO" 3 ). 

In the case of 10 5 Mq IMBHs the constraints are even 
stronger: HoII cannot host 10 IMBHs with m. = 10 s M©, 
unless their efficiency is 77 £ 10 -6 . 

From these results, we derive the following upper limit 
for disc IMBHs: 



n. < io" 5 r2( 



Go- 3 ) 



/ m, 

\ 10 4 M e 



(12) 



This approximate formula is valid for 10 — 10 Mg IMBHs. 
We made some checks for more massive BHs (10 6 — 10 7 Mq) 
and we found that it is still applicable; but we cannot verify 
whether it holds also for lower masses. In equation (12) the 
scaling with tj" 1 is entirely due to the fact that we adopted 
the Bondi-Hoyle formula (equation 8) to calculate the lumi- 
nosity. The scaling with m. _1 is also due to the Bondi-Hoyle 
formula and to the fact that Q., scales with the IMBH mass 
(see equation 6). 

The constraints on halo IMBHs are looser by a factor 
of ~ 10 3 . 1000 IMBHs with mass m. = 10 4 M Q (i.e. fi. = 
10 -2 fi(,) overproduce the number of X-ray sources if their 
efficiency is 77 > 10" 3 . Similarly, 1000 halo IMBHs with mass 
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Figure 4. Cumulative distribution of observed and simulated X-ray sources as a function of the X-ray luminosity Lx- The lightly 
hatched histogram reported in all the panels is the cumulative distribution of observed X-ray sources (Kerp et al. 2002; the luminosity 
of HoII X-l is assumed to be in its brightest state, Lx = 2 X 10 40 erg s — , Dewangan et al. 2004). Open and heavily hatched histograms 
refer to simulations. Left panels: initial conditions; right panels: after 6 Myr. First row (from top to bottom): disc4_10 (open histogram) 
and disc4_100 (heavily hatched); second row: disc5_10 (open); third row: DMM4.100 (open) and DMM4.1000 (heavily hatched); fourth 
row: DMM5.100 (open) and DMM5.1000 (heavily hatched). 



m. = 10 5 Mq exceed the number of X-ray sources in HoII 
for r) > 10" 5 . 

Thus, an upper limit valid for 10 4 - 1O 5 M halo IMBHs 
(distributed according to DMM) is: 

which, a part for the normalization, is the same as equa- 
tion (12). However, this constraint is more affected by sta- 
tistical fluctuations, because of the small fraction of IMBHs 



(~ 3 — 5 per cent) which efficiently accrete with respect to 
the total number. 

This upper limit is similar to the analogous constraint 
that MFR found for 10 4 Mq IMBHs distributed according 
to a DMM profile in the Milky Way (Q. < W~ 2 n b for r/ = 
10" 3 ). In this paper, the same upper limit has been derived 
on the basis of more robust simulations, as MFR used a 
simple Monte Carlo treatment for the gas instead of SPH 
simulations. 

As we said, the efficiency of IMBHs to attract clumps 
of gas particles is important in determining the cumulative 
luminosity distribution. But the dumpiness of the simulated 
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disc is function of the free parameter k. In the simulations 
showed in Fig. 4k = 9 x 10 23 g _2//3 cm 4 s -2 , which corre- 
sponds to a gas temperature T gas ~ 60 K (i.e. the average 
temperature of gas in molecular clouds). What happens if 
we change «? Are our results sensitive to this parameter? 
In Fig. 5 we compare the cumulative luminosity functions, 
after 6 Myr, of the runs disc4_100 (where k = 9 x 10 23 g~ 2/3 
cm 4 s -2 ) and K_disc4_100 (where k = 1.5 X 10 25 g" 2/3 cm 4 
s -2 , corresponding to T gas ~ 1000 K). The two distribu- 
tions are similar. The only significant difference is in the 
high-luminosity tail, because the warmer gas becomes less 
clumpy and its sound speed c s is higher (the average value 
of c s around the IMBHs is (c s ) ~ 3 km s" 1 for disc4_100 and 
(c s ) ~ 6.5 km s _1 for K_disc4_100). Then, k does not signif- 
icantly influence our results, at least for a reasonable range 
of values. Of course, higher k imply larger Jeans masses, 
and then larger scale lengths of the clumps. However, these 
differences do not substantially affect the region close to the 
IMBH, where we sample the density p g . 

In Fig. 5 we can also note the effect of adopting for 
Cs the values derived from the simulation (see Section 4.1). 
While most of the distribution is unchanged, the high- 
luminosity tail is strongly affected by adopting the simu- 
lated c s (upper panel) instead of the expected value (lower 
panel) for an unperturbed cloud (i.e. c s = 0.3 km s _1 ). This 
is due to the fact that the highest density clouds around the 
IMBHs are also the hottest. In the case of the run disc4_100 
the luminosity distribution drops at Lx Si 1 .8 x 10 41 erg s~ 1 , 
when the simulated c s is adopted (top panel); while, if the 
sound speed is assumed to be c s = 0.3 km s _1 , the brightest 
sources have Lx ~ 4.5 x 10 41 erg s" 1 (bottom panel). The 
difference is even larger for the run K_disc4_100, where the 
gas has higher initial temperature (T gas = 1000 K) : the high- 
luminosity tail reaches Lx = 4.5 x 10 41 erg s _1 if c s = 0.3 
km s _1 (bottom panel), and stops at Lx — 4x 10 40 erg s _1 
if we derive c s from the internal energy (top panel). 

As a sanity check of the stability of our results we also 
protracted some runs (disc4_10, K_disc4_100, DMM4_1000 
and HI_DMM5_100) for 10 Myr. We concluded that the sit- 
uation at 10 Myr is not statistically different from the situa- 
tion at 6 Myr, with respect to the inferred accretion rate of 
IMBHs, as most of the density evolution occurs in the first 
5 Myr and the IMBH luminosity distribution remains quite 
stable in the next million years (Fig. 6). 



5 ATOMIC HYDROGEN 

MFR showed that in the Milky Way even IMBHs accreting 
atomic hydrogen can become observable as X-ray sources 
and are crucial to constrain the number of IMBHs. Thus, 
we considered atomic hydrogen also for HoII. Unfortunately, 
we cannot simulate a two-phase medium. Thus, in some runs 
we substituted molecular hydrogen with an exponential disc 
of gas. The total mass of gas as well as the scale length and 
scale height have been chosen according to the observational 
values (see Sections 2 and 3.1). 

To calculate the luminosity of IMBHs accreting atomic 
gas, we adopt the same procedure described for molecular 
hydrogen (Section 4), except for the choice of the radius r g . 
In fact, given the lower density of atomic hydrogen regions, 




1 36 1 Q 37 1 38 1 Q 39 ^ Q 40 , Q 41 , q* 

L x (erg s" 1 ) 



Figure 5. Cumulative distribution of X-ray sources as a function 
of the X-ray luminosity Lx, after 6 Myr, for the runs disc4_100 
(open histogram) and K_disc4_100 (hatched). Top panel: c s de- 
rived from the simulations. Bottom panel: c B = 0.3 km s —1 . 
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Figure 6. Cumulative distribution of X-ray sources as a function 
of the X-ray luminosity Lx for the run disc4_10 at the beginning 
of the simulation (dot-dashed line), after 2 Myr (dashed), after 
6 Myr (solid) and after 10 Myr (dotted). The luminosity rapidly 
increases in the first 6 Myr and then stabilizes. 



we adopt r g = 100 pc, much larger than for molecular hy- 
drogen. 

To avoid contamination with IMBHs accreting molec- 
ular hydrogen, we do not consider those IMBHs that are 
inside the volume which should be occupied by the molecu- 
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Figure 7. Cumulative distribution of X-ray sources as a func- 
tion of the X-ray luminosity Lx, after 6 Myr, for the runs 
HI.DMM5.100 (heavily hatched histogram) and DMM5.100 
(lightly hatched). The open histogram shows the observations 
(Kerp et al. 2002; Dewangan et al. 2004). 

lar hydrogen disc (i.e. within a scale length R mo i = 500 pc 
and a scale height z mo i— 60 pc). 

Fig. 7 shows the results for the run HI_DMM5_100, i.e. 
the case where the contribution of atomic hydrogen has been 
found to be more important. It is apparent that IMBHs 
accreting atomic hydrogen (heavily hatched histogram) have 
luminosity Lx < 10 37 erg s . Thus, these sources are below 
the observational threshold of Kerp et al. (2002), and are not 
relevant for our purposes. 

This is not in contradiction with the results of MFR. 
In fact, also the luminosity distribution of IMBHs accreting 
atomic hydrogen in the Milky Way drops at Lx = 10 37 erg 
s _1 . The difference between MFR and this paper consists 
merely in the data, which in the case of the Milky Way 
are complete down to 10 36 erg s _1 , rather than to 10 37 erg 
s _1 as in HoII. In addition, most of Galactic sources have 
been identified as low mass X-ray binaries, high mass X-ray 
binaries and other classes of known sources, excluding the 
possibility of being IMBHs. 



6 SUMMARY 

In this paper we derived upper limits on the density of 
IMBHs in HoII, by assuming that they can accrete dense 
gas and emit as X-ray sources. In particular, we ran SPH 
N-body simulations of a disc dwarf galaxy having the main 
observational properties of HoII and hosting a halo or disc 
population of IMBHs. We compared the luminosity distri- 
bution derived by these simulations with the observational 
data by Kerp et al. (2002). A similar technique has been 
previously applied to the Milky Way by MFR. 

The target galaxy of the current paper, HoII, was chosen 
for its rich content in gas and X-ray sources, as well as for its 



small mass (if compared to the Milky Way), which enables 
us to reach very high mass and spatial resolution if compared 
to MFR. The improvements with respect to MFR are not 
only in the resolution, but also in the treatment of the gas 
(this paper is based on SPH simulations, whereas MFR used 
a Monte Carlo approach) and in the presence of a 'living' 
halo (MFR adopted a rigid halo). 

MFR assumed that IMBHs, if exist, are a (more or less 
concentrated) halo population. This might be the case if 
IMBHs are the remnants of massive population III stars 
formed in minihalos. But other formation mechanisms, such 
as the runaway collapse in young clusters (Portegies Zwart 
& McMillan 2002), suggest that IMBHs can also be a disc 
population. IMBHs could reside in discs even if they are 
remnants of population III stars formed at relatively low 
redshift (Krolik 2004; Jimenez & Haiman 2006). Thus, in 
this paper we consider both the case where IMBHs are a halo 
population (following a DMM profile) and the case where 
IMBHs are a concentrated disc population. 

If IMBHs are a disc population, we find an upper limit 
SI. < 10~ 5 f2;, 77^ 1 1 (where 773 = 77 /10~ 3 and nn — 
?n./10 4 M Q , equation 12). 

In the case IMBHs follow a DMM profile, the upper 
limit can be written as S7. < 10 _2 f2t, mj 1 (equation 13), 
similar to the upper limit found by MFR for the Milky Way. 

Thus, limits on disc IMBHs are a factor of 1000 stronger 
than limits on halo IMBHs. This implies that only 1 IMBH 
with mass 10 4 M Q can be hosted in the disc of HoII, unless 
the radiative efficiency is lower than 10~ 3 . Such limits have 
important consequences on models of IMBH formation such 
as the runaway collapse scenario. 

In summary, the results presented here confirm the re- 
sults by MFR, and strengthen the upper limits derived from 
the Milky Way. Unfortunately, this method suffers of a num- 
ber of uncertainties. First of all, the radiative efficiency r\ 
was left as a free parameter (even if we adopted r\ = 10 -3 
as a reference value for the figures). In addition, our method 
does not account for the effective duty cycle of such IMBHs. 
As suggested by Krolik (2004), the duration of the high- 
luminosity state of these sources must be very short (less 
than 1 Myr), because the IMBH rapidly heats the surround- 
ing gas, quenching the accretion rate. This effect cannot be 
considered in our simulations, and should be further inves- 
tigated. 
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